home *** CD-ROM | disk | FTP | other *** search
/ SGI Hot Mix 17 / Hot Mix 17.iso / HM17_SGI / research / lib / obsolete / gauss.pro < prev    next >
Text File  |  1997-07-08  |  954b  |  56 lines
  1. ; $Id: gauss.pro,v 1.2 1997/01/15 04:02:19 ali Exp $
  2. ;
  3. ;  Copyright (c) 1991-1997, Research Systems Inc.  All rights
  4. ;  reserved. Unauthorized reproduction prohibited.
  5.  
  6.  
  7. function gauss, a ;
  8. ;+
  9. ; NAME: 
  10. ;    GAUSS 
  11. ;
  12. ; PURPOSE: 
  13. ;    Gauss returns the cutoff value v such that
  14. ;
  15. ;        Probability(X > v) = a,
  16. ;
  17. ;    where X is a standard gaussian random variable.
  18. ; 
  19. ; CATEGORY:
  20. ;    Statistics.
  21. ;
  22. ; CALLING SEQUENCE: 
  23. ;    Result = GAUSS(A)
  24. ;
  25. ; INPUT:
  26. ;    A:    The probability for which a cutoff is desired.
  27. ;
  28. ; OUTPUT: 
  29. ;    The cutoff value if a is beween 0 and 1 inclusively. Otherwise, -1.
  30. ;-
  31.  
  32. if (a gt 1 or a lt 0) then return,-1  
  33. if a eq 0 THEN return, 1.e12
  34. if a eq 1 THEN return,-1.e12
  35.  
  36. if (a gt .5) THEN BEGIN
  37.   a = 1-a
  38.   adjust = 1
  39. ENDIF ELSE adjust = 0 
  40.  
  41. BELOW = 0
  42. UP    = 1.0
  43.  
  44. while gaussint(UP) LT 1.0 - a DO BEGIN
  45.  Below = UP
  46.  UP = 2*UP
  47. ENDWHILE
  48.  
  49. x = pd_bisection([1.0-a],'gaussint',Up,Below)
  50. if adjust THEN BEGIN
  51.   a = 1 - a
  52.   return, -x
  53. ENDIF ELSE return,x
  54. end
  55.  
  56.